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_. ' Abstract 

« ; 

^ I We investigate some classical evolution model in the discrete 2 + 1 space-time. A map, 

giving an one-step time evolution, may be derived as the compatibility condition for 

some systems of linear equations for a set of auxiliary linear variables. Dynamical 

pg ■ variables for the evolution model are the coefficients of these systems of linear equa- 

^ , tions. Determinant of any system of linear equations is a polynomial of two numerical 

C^^ ' quasimomenta of the auxiliary linear variables. For one, this determinant is the gen- 

^^ I crating functions of all integrals of motion for the evolution, and on the other hand it 

^T^ • defines a high genus algebraic curve. The dependence of the dynamical variables on 

^—v , the space-time point (exact solution) may be expressed in terms of theta functions on 

0^\ ' the jacobian of this curve. This is the main result of our paper. 

o^: 

-*— > ■ 

Ch ■ 1 Introduction 

J^ [ In this paper we give an exact solution of a classical evolution model in discrete 2-1-1 

O ■ space-time. This model was formulated in ||l|. The map of dynamical variables, governing 

• • . the one-step evolution, was derived as the compatibility condition for two sets of linear 

.^^ I relations, associated with the usual graphical representation of left and right hand sides of 

/\ ' the famous Yang - Baxter equation. Such form of the zero curvature condition generalizes 

j^ ■ the usual concept of the local Yang - Baxter equation as the zero curvature condition for 

3d lattice models. Linear variables may be assigned either to the vertices or to the sites 

of the two dimensional auxiliary graphs (these two assignments are the dual ones), and 

the coefficients of the linear relations - nothing but the dynamical variables - are to be 

assigned to the vertices. Main feature of the map of the dynamical variables is that it is 

canonical with respect to local Poisson brackets and thus may be easily quantized |y], so 

that the auxiliary linear systems exist even in the quantum case |Q]. 

Evolution models arise when one considers flat regular graphs, formed by straight lines, 
on a two dimensional torus. The only demand is that the graph must contain Yang-Baxter 
triangles, so that a simultaneous "bypass of some lines through appropriate vertices" will 
restore the geometry of the graph. Thus, with the map associated with the Yang - Baxter 
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type triangles, there appears the map of the dynamical variables, assigned to the vertices 
of the primary lattice, into the set of the dynamical variables, assigned in fact to the same 
lattice. This gives the one-step evolution of the discrete 2 + 1 evolution model. 

Simple square lattice does not fit the demands, because of it does not contain the 
triangles. The simplest two dimensional graph with the properties demanded is the so 
called kagome lattice Q (see the figures below). This is not something strange, the 2d 
kagome lattice is nothing but a section of a regular 3d cubic lattice by an inclined plane. 

The governing map was derived for the open linear system (i.e. the number of the 
equations is less then the number of auxiliary linear variables), but nevertheless the linear 
systems may be written for the 2d lattice with the toroidal boundary conditions. Also, 
in general, dealing with the linear equations, one may impose on the linear variables 
quasiperiodical boundary conditions in both directions of the torus with two numerical 
quasimomenta. In this case the number of the auxiliary linear variables coincides with the 
number of the linear equations and one may ask for the admissibility of the linear system, 
i.e. for the zero value of corresponding determinant. Because of the evolution arises as a 
simple compatibility of two similar linear systems, the admissibility of the primary system 
provides the admissibility of the evaluated one. Therefore the determinant is at least 
an ideal of the evolution. Moreover, being normalized appropriately, formal determinant 
J{A, B) of the linear system as a polynomial of the quasimomenta A and B is conserved by 
one-step evolution map, i.e. is the generating function of the integrals of motion. All these 
remains valid and in the quantum case, where J{A^ B) is an operator-valued functional. 

In our classical case and finite spatial size of the two dimensional lattice, equation 
J{A,B) = gives a finite genus g algebraic curve T, so that the integrals of motion are 
interpreted as moduli of T. With F and a bit of additional information concerning the 
initial state given, the system of the auxiliary linear variables may be easily parametrized 
as the meromorphic functions on T in terms of the theta functions on Jac T. Doing this, 
we get at once the parametrization of the dynamical variables in terms of theta functions 
on Jac r and obtain the exact solution. 

Perhaps it would be expedient to discuss several 3d discrete integrable models from the 
point of view of their linear systems and indicate the place of the model being considered 
among them. Spatial nature of any 3d integrable model means that geometrically linear 
variables are assigned to several elements of 3d cubic lattice. This assignment gives a 
type of the linear system. Consider three main scenarios corresponding to three relative 
integrable models. 

The first, most simple type [g, y, Q, §]: let the vertices of the cubic lattice Z^ have the 
coordinates p = (a, b, c), a,b,c G Z. Consider the even sublattice of it, Z^^^^^i p = (a, b, c), 
a + b + c = even. Points Z^^^^ triangulize the three dimensional Euclidean space into 
the following convex bodies: set of octahedra and sets of two types of tetrahedra (up to 
regular translations). Assign the auxiliary linear variables to the vertices of Z^^^^. Linear 
equations are assigned to the triangles of the triangulation described. Primary set of the 
linear equations maybe chosen for the triangles - sides of one of the tetrahedra. System 
of coefficients of the linear equations gives the tau function for Hirota's discrete bilinear 
equation (the ^^ octahedron equatiorH^ from this naive geometrical point of view). 

The second type: assign the auxiliary linear variables to the facets of Z^. Linear 



^ "Kagome" is not a name, it is a kind of Japanese mats. 
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Figure 1. Labelling of the kagome lattice. 



relation are to be written for each edge surrounded by four facets. This corresponds to 
Korepanov's block-matrix models M. Hirota and Hirota-Miwa' equations are the simplest 



12 



compatibility conditions for Korepanov's linear problem [11, |1C 

The third type: dealing with the cubic lattice, most obvious scenario is to assign the 
linear variables to all the vertices of it (or, correspondingly, to the sites of the dual lattice). 
Each linear equation corresponds to a square (any face of the elementary cube). This is our 
case. Independent Lagrangian - type variables are a triplet of "tau - functions" , obeying 
the system of the "cu6e equations" . Taking a section of the cubic lattice by an inclined 
plane, we at first get the kagome lattice geometrically, and the linear system for it as a 
part of whole linear relations secondly. Details may be found in |j^, |2|, |3| . The advantages 
of this approach (Poisson structure, quantization etc.) were mentioned at the beginning 
of this introduction. 

Few remarks concerning the section by the inclined plane and the evolution. For the 
cubic lattice p = (a, b, c), a,b,c G Z, the sections mentioned are the planes a + b + c = t = 
const. Equations of motion normally may be solved so that all the dynamical variables 
for a + b + c = t + l are expressed from the dynamical variables for a + b + c = t. This 
gives the natural notion of the evolution as the map from t to t + 1. As it was mentioned, 
geometrically the section a + b + c = t is the kagome lattice. 

This paper is organized as follows. First, we recall the formulation of the model, 
describing the dynamical system and defining the evolution. Second, introducing the 
linear system, we define the generating function for the integrals of motion. All these are 
based on the results of [||, ||, and we rather enumerate the facts. In the fourth section we 
analyse the curve, parametrize the auxiliary linear variables and derive the expressions for 
the dynamical variables. In the final section we discuss possible applications of the results 
obtained. 
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Figure 2. (a, b) - triangle. 

2 Discrete evolution 

System of the dynamical variables assigned to a finite M xM kagome lattice (see Figure [l|) 
is a set of 3 M^ pairs 



[yij,a,b , v\'j,a,b] , j = 1,2,3, a,b G Zm ■ 



(2.1) 



Arrangement of the triangles of the kagome lattice is shown in Figure ||. Consideration 
of the lattice on a torus implies the periodical boundary conditions for the dynamical 
variables 



^j,a+M,b — ^j,a,b+M — ^j,a,b , ^j,a+M,b — '^j,a,b+M — ^j,a,b ■ 



(2.2) 



Geometrically the indices (j, a, b) are assigned to the vertices of the kagome lattice, so 
that for a,b given {j,a,b), j = 1,2,3 mark three vertices of a definite triangle, as it is 
shown in Figure 0. Whole kagome lattice may be obtained as the up-down and left-right 
translations of the triangle (a, b). It is supposed that a increases to up and b increases to 
right. 

Impose the following Poisson structure on the set of the dynamical variables: 



{ ^j,a,b , '^j,a,b } — ^j,a,b '^j,a,b , 



(2.3) 



and any other Poisson bracket is zero. Remarkable is the locality of the dynamical vari- 
ables. 

Evolution of the system is governed by the fundamental map R. Consider one isolated 
triangle (a, 6). Define map 



R : [uj,Wj] 1-^ [u',w'] j = 1,2,3, (a, 6) implied to be fixed. 



(2.4) 
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Figure 3. Geometrical representation of the map R. 



by the following relations: 



Wi 



K2 Ui U2 W2 



Ki U^ W2 + K3 U2 W3 + K^ Kg U^ Wg 

W^ W2 + U3 W2 + K3 U3 W3 

Wo 



(ii) < 



u'. 



W9 



Ul U2 U3 



U2 U3 + U2 W-j^ + K^ U^ W-j^ 

W^ W2 W3 
Wi W9 + Uo W, + Kn Uo W 



J3 VV2 



3 "3 ™3 



(2.5) 



[m) 



Wo 



Uo Uo + Un Wi + Ki Ul W 



■'2 U3 



1 Ul w^ 



K2 U2 W2 W3 



K^ Ul W2 + K3 U2 W3 + K]^ K3 U^ W3 

K]^ 2 3 ^^6 arbitrary numbers (not the dynamical variables). 

Proposition 1. The map R, l \2.3^ , conserves the local Poisson structure, 

{Uj,\Nj,} = 6jj,Uj\Nj 4^ {u'j,\n'j,} = Sjjiu'jVj'j, 

i.e. R is the canonical map. 



(2.6) 



Geometrically R may be interpreted as the map from one Yang - Baxter triangle to 
another, as it is shown in Figure ^, i.e. as a bypass of a line through the opposite vertex. 

Turn now to the whole lattice. We will distinguish three types of lines with respect to 
their slopes by the letters x, y, z as it is shown in the Figure Q and enumerate the lines, 
Xa,yi3,z^, CK)/3,7 G Zm- "Spectral parameters" Kj^a,b actually depend on numbers of the 
corresponding lines. 
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Figure 4. Accessories of the lattice: lines. 



Evolution of the lattice is the simultaneous shift of all lines of one type in one direc- 
tion. The result of such shift is the application of R to each (a, b) triangle and some 
re- numeration of the images of u'- ab^^'j ab- This re- numeration depends on what type of 
vertices we leave immovable. Choose the vertices of type {l,a,b) motionless, i.e. we shift 
all X lines to the north-east direction. Define the functional map U, acting on the space 
of functions of the dynamical variables [uj^a,b,'^j,a,b]'- 



(U o /) ( Uj^a,b , '^j,a,b ) = /( U* o Uj^a,b , U* o w 



j,a,b 



(2.7) 



where 



U* o ui „ . = u'l „ f, , U* o wi , 



w 



l,a,fe 5 



U* o U^,a,b = U2,a+1,6 , U* o w^^^f, = Vj'2^a+l,b > 



(2.^ 



U* o u 



3.a,b 



U3,a,b+1 , U* o Wg^^^, 



W 



3,a,fe+l ! 



where, for example, u'^ ^^ means that we take the expression for u'^^ from ( |2.5D and add the 
indices a, b to each Uj, Wj there. Indices (a + 1, b) and (a, b + 1) in the second and third 
lines of (^]^) are the re-enumeration mentioned above. 

Obviously, U conserves the Poisson brackets, and thus it is the canonical map. U is 
identified with the one step discrete evolution, so that if 



/( Uj,a,b , Wj,a,b ) = f{to) , 



then 



/(io, t) = (U*o/)(to). 



(2.9) 
(2.10) 



3 Linear system and the integrals of motion 

Map R ( p. 5]) was "derived" in Q as a zero curvature condition for a system of linear 
equations. There are at least two ways to define the linear system, and here we use the 
co-current form according to the terminology of [|l|, §]• 
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fd R 




Figure 5. Linear variables of the map R. 

We start from the linear system for isolated Yang - Baxter triangles and explain how 
the map R appears. First, in addition to the vertex variables Uj,\nj,kj introduce eight 
auxiliary variables ipa-.-ifh living in the sites of the 2d graphs as it is shown in Figure ^ 

Please note that the variables ipi,...^pg are the same in the left and right hand sides of 
Figure g and belong to the equivalent open cells. Left and right hand side graphs differ 
by iph and (fa- Consider now two sets of linear relations: for the left hand side graph 







def 
Jl = "Pc 



fe^l + VfeWi + (^rfKiUiWi , 







n d,ef 

h = "Ph 



V>d^2 + Vb'^2 + ^fl^2^2'^2 > 

and for the right hand side graph 

fa u'l + fb w'l + ff 1^1 u'l w'l , 
= /2 =' fc - 'Pe^'2 + fg W2 + Va 1^2 ^2 ^2 ' 
= /s = ^a - V^e U3 + V5/ W3 + C^^ K3 U3 W3 . 



(3.1) 



= A =' f. 



h = 'Pc 



(3.2) 



Each linear expression fj and /' is assigned to j'-th vertex of the left and right graphs 
respectively. Into each fj and /' there involved four linear variables from the sites sur- 
rounding the vertex. Coefficients in fj and /^ are the vertex variables, assigned with j-th 
vertex. Arrangement of the vertex variables and the site variables with respect to the 
arrows of the lines is the same for any vertex (just look attentively at the linear relations 
and the figure). Each expression fj and /' becomes the equation, fj = 0, /' = 0, and thus 
two sets of linear equations appear. 



Theorem 1. Two systems, (3.1) and (3.<) are linearly equivalent (after excluding extra 
linear variables iph and ipa) ijff u',w'- are connected with Uj,\Nj via ^A 
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'P3,a,b 



'P2,a-l,b 

Figure 6. Accessories of the lattice: cocurents. 



Turn now to the kagome lattice on the torus. Linear forms fj,a,b and and corresponding 
linear equations fj,a,b = are to be introduced for all vertices of the lattice. 







„ def 

Jl,a,b = ^3,a+l,b 

+ V3,aJ 



r def 

J2,a,b = Vl,a,6 



- 9^2,0,6 • ^l,a,b + fl,a,b ' '^l,a,i 
-1 ■ 1^1 Ul,a,6Wi afc , 



^3,a,b+l • U2,a,6 + 9^3,0,6 ' 'J^2,a,b 
+ ^2,a-l,b ■ K2U2nb^2,a,b , 



(3.3) 



def 



= f3,a,b = ^3,a+l,b — ^l,a,b ' yi3,a,b + V^2,a,fe-1 ' '<'^3,a,b 
+ 'P3,a,b ■ K3 U^^a,b '^3,a,b ■ 

What these writings mean. The linear objects ^jab^ appeared in these relations, are 
assigned to the sites of the kagome lattice as it is shown in Figure |6|. Each linear expression 
fj,a,b is assigned to {j, a, 6)-th vertex of the kagome lattice, and into fj,a,b there involved 
four linear variables from the sites surrounding the vertex. Coefficients in /j,a,fe are the 
vertex variables, assigned with (j, a, 6)-th vertex. Arrangement of the vertex variables and 
the site variables with respect to the arrows of the lines is the same for any vertex. Each 
expression fj,a,b becomes the equation, /j,a,fe = 0, and thus the set of linear equations 
appears. 

The linearity of cp- ^ ^ allows one to impose the quasiperiodical boundary conditions for 
them: 



<fa,a-M,b 



<^a,a,b ■ A ■, ^a,a,b-M — <^a,a,b ' B 



(3.4) 



The coefficients of the linear system form a 3Af^ x 3M^ matrix L, depending on the 
dynamical variables and the quasimomenta A, B. The following set of propositions was 
proved in Q |2[ : 



Proposition 2. The determinant of Ij is a Laurent polynomial of A,B, 
detL = J2 J^^P^"" •^^' 



(3.5) 
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where domain U is described below. Being normalized in any way, 

J{A,B) = p^, (3.6) 

functional J{A, B) is the generating functions for the integrals of motion, 

J{A,B) = J2 JaA^^^)^'" •^^ (3-7) 

a,/3en 

{U o J„_^)(u,w) = Ja,p{u,w) . (3.8) 

Proposition 3. Domain 11 in the decomposition of J{A, B) is the following hexagon: 

n : -M < a < M , -M < /? < M , -M < a + (3 < M, (3.9) 

where M is the spatial size of the kagome lattice. According to the Riemann-Hurwitz 
theorem, the genus of the curve T : J{A, B) = is 

g = 3M^-3M + l. (3.10) 

Complete number of the integrals of motion is 3M^ + 1, and one can choose exactly SM"^ 
involutive between them. 

Perimeter of the hexagon 11 is formed by 6M points, 

JM-n,n ) J-n,n-M 5 Jhl-n ; J-M,M-n ; Jn-M,M ; Jn,-M ; (^-H) 

where n = 0, ...,M. These perimeter integrals are not independent. Let 

^/3 = (-)^^ n ^i.A-^ia' (3.12) 



^7 = n ™i 



(T,7 ™2,(T,7 • 



Each of these expressions corresponds naturally to its line, X^ to Xa etc., see Figure Q. 
The perimeter integrals ( 3.11| ) are some the symmetrical polynomials of Xa, or Y^, or Z^. 



The Poisson brackets for Xa,Yp, Zy Va, (3, 7 are 

{^0,^/3} = XaYp , {Yf3,Zy} = Yf^Z-f , {Zy,Xa} = ZjXa. (3.13) 

It is useful to extract common convolutive parts from Xa,Y/^,Zy, 

X^ = X ■ j{Xa) , Yp = Y ■ j{Yp) , Zy = Z ■ j{Zy) , (3.14) 

where 

{X,Y] = XY , etc. (3.15) 
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Now we can rewrite J{A, B) in its final invariant form: 

B\'^ 
1 



a,/3,7en' 



(3.16) 



where H' is three squares of the cube: 

n' : < a,/3,7 < M , at least one of a,/?, 7 iz zero 



(3.17) 



Now we may describe the complete set of SM"^ involutive integrals: g = SM"^ — 3M + 1 
functionals ja,f3,-y, corresponding to the inner points of XT', 3M — 3 independent projective 
j{Xa), j{Yfs) and j{Z.y), one X -Y ■ Z (this is the center for X,Y, Z), and any finally - 
any other single function of X, Y, Z. 

Running ahead, for what purpose else one needs the perimeter integrals ( 3.11 )? 



J{A,B) = is an algebraic curve T, and we will be interested in the divisors (A) and 
(-B) of the algebraic functions A = A{P) and B = B{P), P G F. In general the perimeter 
integrals describes the divisors {A) and (S), and after a bit cumbersome calculations we 
have obtained the following description: 



and 



(^)c 



(^)c 



A 
A 
A 

A 



A 



Xa , a ^ Zm 



, B = Yp, (3£Zm, 

A Xr, 



°"' B 



(^^2^.3)^ 



, a G Zm , 



00 



B 



— ) Y(s, PgZm 



(3.18) 



iB\ 



iB)c 



B 



B 



B 



B 



, — = Xa , a £ Zm 



0, A 



A 

°"' B 



1 



{H1H2Y' Z^ 



7 e Zm 






(K2K3) 



aGZ, 



M 



00 , A = —- 



7 e Zm , 



(3.19) 



In these formulae it is supposed that the spectral parameters Hj are the same for all lattice, 
but the structure of Xa, Yp and Z^, associated with the lines Xa, yp and z^ makes obvious 
the situation when Kj depend on the lines numbers. 
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4 The curve and solution 

Turn at last to the algebraic geometry. We deal with the algebraic curve F, defined by 
J{A, B) = {). A, B and C = A/B are the meromorphic functions on F, their divisors 
are already obtained in the previous section, and ( 3.18 , pyi9 ) may be rewritten in decent 
notations as 

(^) = E iPl - p^J + E (n^i - Pi, 

a p 

(4.1) 

a 7 

where the notion of the points P^^, P^, and P^ comes from ( p.l8 , 3ri9| ) transparently. 
Remarkable feature of these divisors is the natural correspondence between P and the 
lines Xa, yf3, z^- 

Return now to linear system (|3.3| ). Solving it, one may put one of ^j^a,b to be unity, then 
all other linear variables become some meromorphic functions on T. Our next interest is a 
common pole divisor of all ^j,a,b- Denote this divisor as V, its degree may be calculated: 

deg:D = 3M2 + 1 . (4.2) 

This calculation is based directly on the dimension of linear system and needs no com- 
ments. The set ^Pj,a,b is the unique solution of the linear system, so due to the Riemann- 
Roch theorem the dimension of the linear space of the holomorphic functions with the 
pole divisor V is 

dim LCD) = degV - g + I = 3M + 1 . (4.3) 

This means that one may restore ^j^a,b as a function on T via just 3M points of positive 
^j,a,b^ 

{ipj,a,b) +V = vfl, + V'^^,^,, degVfl, = g, deg P^- ,,, = 3 M , (4.4) 

where positive T)- j ^ may be restored unambiguously. Obviously, the meaning of the 
quasimomenta A and B implies that we may choose the lines Xa^ , yp^ and z^,, , where the 
quasimomenta appear, in any way. Therefore Pp, p = {j,a,b), are all governed by the 
same points as form the divisors (A) and {B). 
The same is valid for V: 

f) = p(0) + f)' ^ (45) 

where P', deg P' = 3M is also governed by the points of {A) and {B). Thus for any 
p = (j, a, b) the decomposition arises: 

(<^p) = ((^p)o - ((^p)oo = V^p^ - PW + P'p - V' , (4.6) 

such that we can trace the following simple part of (v^p): 

Pp = P'p - V' . (4.7) 
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P" : ^c 





p+ : Pp+ = Pp_+P+-P- 



Figure 7. Sites separated by a segment of line x between two vertices and theirs divisors. 



Before we give concrete expressions for cpp recall a couple of notations of the algebraic 
geometry, see for example [|l^]. For a given curve F with normalized holomorphic one- 
forms u! it is defined its jacobian Jac F and the theta - functions on it. We'll use the 
conventional notations for F i— > Jac F: for D, deg T> = 0, let 



liV) 



u 



Also for P, Q E F denote their prime form as E{P, Q), 



E{p,Q) ~ e^ 



LO 



(4.8) 



(4.9) 



where the subscript 5 means the nonsingular odd theta characteristic of the curve (see 
|T3|). Let the same symbol E stands for a product of the prime forms 



E{P,V) = H E{P,Q), V > 0. 
Qev 



(4.10) 



Theorem 2. With the notations introduced, the expression for ipp as the meromorphic 
function on F with the divisor given by ^dX,4--^i) is 



,p, _ (0) e(z + /(p-Po + Pp)) i?(p,(Pp)o) 

^9 - MP) - ^P e(z + /(P-Po)) E{P,{VpU[ 



(4.11) 



where P G F, vector z G C^ and Pq € F are some auxiliary parameters in the parametriza- 
tion ( 4-11 ), ^p are constants (i.e. do not depend neither on P nor on z), 0(z + /(P — Pq)) 
corresponds to g poles from T>^^', and so on. 



Actually eq. ([4.11 ) is the consequence of the Riemann-Roch theorem, see ||9|, p^ . 

Further it is simpler to cancel G(z+/(P— Pq)) from all yjp and deal with the holomorphic 
with respect to z functions. 

Describe now the explicit way of assigning the divisors Pp ( [4.7| ) to the sites p = {j, a, h) 
of the kagome lattice. 

To each line Xa, yp and z^ of the lattice the pair P+ — P7 is assigned, see ( [l.ip and the 
remarks at the beginning of this section. Let a segment of oriented line x separates two 
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X : V, 




X : D, 
Figure 8. Divisor notations around a vertex. 

sites, p~ on the left and p+ on the right according to the orientation of x, see Figure ^. 
According to our previous considerations, the divisors Pp+ and Pp- , assigned to p"*" and 
p~ respectively, obey 

where P^ and P~ is the pair assigned to line x. Starting from any site pg on the lattice 
and using this procedure, we may define all Pp up to the divisors of the algebraic functions 
[A) and {B). Such system of the divisors was introduces by I. Korepanov in |^], i.e. the 
divisor rules we've obtained (Figure 0, eq. ( [4.12| )) coincide formally with that of the second 
type linear problem (see the introduction). 

Introduce also a divisor of the edge of x separating p"*" and p~: 

V, = Pp_ + P+ - Po = V + P^ - Po , (4.13) 

where Pq is the same point as in ( [4.11| ). The meaning of T)^ is following: 

(0) 

^ (z = 5 - /(PJ) = \ . (4.14) 

Consider further a vertex formed by two lines x and y surrounded by the sites a,b,c,d 
as it is shown in Figure S. The linear relation for it is 

(fa -(/?&• U + V^c ■ W + (y9d • KUW = . (4.15) 



With the parametrization ( 4.11 ) this relation must be the identity both in P and z, so 
that neither u nor w depend on P, but u = u(z), w = w(z). This allows one to find u(z) 
and w(z) immediately. 

Eight divisors correspond to eight elements of the vertex. For the site divisors we use 
obvious notations Va, V^, Vc and P^. Two lines involved, x and y, are divided by the 
intersection point into four edges with the divisors: V^ and V^i for x line, and T>y and 
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Vy' for y line. According to divisor rules ( ^.121) , the relations for all these eight divisors 
may be written in the following form: 



^a = 


= ^. 


- p: + ^0 = 


= -^y 


-n 


+ ^0, 


T^b -- 


= ^.' 


- P^ + ^0 = 


- ^y 


-Py 


+ ^0, 


-De -- 


= ^x 


- P^ + ^0 = 


- ^y 


-Py^ 


+ Po, 


-Dd -- 


= ^.' 


- P^ + ^0 = 


- ^y 


-Py 


+ Po- 



(4.16) 



Testing ( |4.15|) for the quasiperiodicity on z and taking zeros and poles of u(z) and w(z) 
(relation (|4.14D is useful for it), we get unambiguously 

B(z + l{Vx>)) B(z + l[Vy)) 

Extra constant parameters obey 



y,W J^(P+,P+) ^f^ EiP+,p- 



^W ' EiP.-,P+)' ^(0) ^ E{P.-,Pi)' 

^_ EiP,+ ,P+) ^__ E{P-,P+) 

^1°) "'' EiP^\Py~)' ^(0) "'°-i?(p-,p-) 



(4.18) 



These constants of the normalization, uq and wq, are to be used to parametrize the "gauge" 
integrals of motion X^, Yj^, Z^, and thus have no dynamical sense. 

Actually ( [4.15 ), as an identity on z, is a combination of two Fay's identities (see |jl3|). 
In the parametrization of u and w the theta functions are assigned naturally to the edges. 
These thetas are nothing but the triplet of tau- functions and thus solve a system of trilinear 
relations, see |^]. 

With the parametrization in terms of the divisors, the evolution becomes simple. It 
corresponds to the geometrical shift of the lines generating the change of the divisors of 
the sites. In general the lattice is heterogeneous and non-equidistant and thus there the 
formula for u,w is rather geometrical, because of "time" means a concrete geometrical 
configuration. But M-step evolution S = U has the common form for all u,w: 

[u,w] = [u(z),w(z)] A [u(z + T),w(z + T)], (4.19) 

where the vector of the periods 

T = I{T) , (4.20) 

is the same for all three "time" directions: 

/ J 2^CK ^oc / J Vol Vol / j Zql -Sq ' \ ' / 

Q O Q 

where, surely, the equahty signs mean the equivalence of the divisors. 
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5 Conclusion 

In this conclusion we would like to discuss possible applications of all the considerations 
made in this paper. Usually the algebraic geometry gives rather formal results, which are 
hard to apply to calculate something "physical". Nevertheless one may use the consid- 
erations above in order to make some important conclusions for the quantum evolution 
model @. 

The quantum evolution model is based on the quantum analogue of the map R ( |2.5D 
for the local Weyl algebra replacing the local Poisson algebra 



Uj ■ Wj = q Wj ■ Uj , (5.1) 

where q is the commonly accepted parameter of the quantum deformation. The form of 
the map R and the linear system in the quantum case differ very slightly form ( ^.5| ) and 
(3.1,3.2), see [0] for the details. The determinant in the quantum case is well defined and 



also gives the complete set of the quantum integrals of motion. With a great effort we do 
not fall into speculations concerning a quantum jacobian. 

Instead, turn to the case when q is A^-th root of unity. For these q the Weyl algebra (|5.l[) 
has the finite dimensional representations, so that A^-th powers u^ and w^ are centres, 
i.e. the parameters of the finite dimensional representation. The map for u^ and w^ 
coincides with the functional map ( p.Sp up to kj i— > k^ , i.e. exactly the case considered. 
Actually the quantum map factorizes into a finite dimensional part (which is the vertex 
R matrix for a Zamolodchikov - Bazhanov - Baxter - type model, giving the Boltzmann 
weights for a statistical mechanics modelling) times the functional part for A^-th powers. 
The evolution operator factorizes also in this way. The quantum S matrix is a high power 
of the quantum one-step evolution operator, and so in order to factorize it one has to carry 
out all the functional parts outside. Doing this, one changes the parameters of the finite 
dimensional one-step evolution operators. This situation corresponds to the consistent 
changes of the Boltzmann weights parameters from one time layer to the other. In two 
dimensions the simplest such models are known as the checkerboard models, and the 
simplest three dimensional chess model was described in [^, 15 1. Note, the parameters 



u^^,\N^,K live on the lattice, i. e. Uy, Wy, Ky are assigned to the vertex V as well 
as the operators uy^'^v- Actually, for given spatial size of the kagome lattice M, one 
may consider another effective spatial size M' of the lattice of parameters, such that M 
is divisible by M'. Small M' correspond to the small heterogeneity of the lattice, and 
M' = 1 corresponds to the spatial homogeneity of the parameters. 

In general the results of this paper would help one to parametrize the consistent evolu- 
tion of the parameters of the finite dimensional one-step evolution operators or the transfer 
matrices. 
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